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Abstract 

The variation of total solar irradiance (TSI) has been measured since 1978 
and that of the spectral irradiance for an even shorter amount of time. Semi- 
empirical models are now available that reproduce over 80% of the measured 
irradiance variations. An extension of these models into the more distant 
past is needed in order to serve as input to climate simulations. Here we 
review our most recent efforts to model solar total and spectral irradiance on 
time scales from days to centuries and even longer. Solar spectral irradiance 
has been reconstructed since 1947. Reconstruction of solar total irradiance 
goes back to 1610 and suggests a value of about 1-1.5 Wm~^ for the increase 
in the cycle-averaged TSI since the end of the Maunder minimum, which is 
significantly lower than previously assumed but agrees with other modern 
models. First steps have also been made towards reconstructions of solar 
total and spectral irradiance on time scales of millennia. 
Key words: solar irradiance, solar magnetic fields, solar-terrestrial 
relations, solar variability 
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1. Introduction 



Solar irradiance is tlie total energy flux (or energy received per unit 
area and time) at the top of the Earth's atmosphere. It is strongly wave- 
length dependent: less than 8% of the solar energy is emitted at wavelengths 
below 400 nm, more than 60% come from the wavelengths between 400 



and 1000 nm and rough 



Solanki and Unruh 



1998 



y another 30% from the longer wavelengths (e.g., 



Krivova et al 



20061) 



The small contribution of the UV part of the spectrum is, however, com- 
pensated by the spectral dependence of the transmission of the Earth's at- 
mosphere. Thus solar radiationbelow 300 n m is almost completely absorbed 



in the atmosphere (see, e.g.. 



Haighl . 120071 . and references therein) and is 



important for the chemistry of the stratosphere and overlying layers. In par- 
ticular, radiation in the Ly-ct line (121.6 nm) and in the oxygen continuum 
and bands between 180 and 240 nm controls production and d e struc t ion o 



ozone (e.g. 



Frederickl . 



1977 



Brasseur and Simon 



1981 



Haighl. 



1994, 



2007 



Fleming et al 



1995 



Egorova et al. 



2004 



Langematz et al. 



20051 ). Solar 



UV radiatio n at 200-350 nm is the main heat source in the stratosphere and 



mesosphere ( iHaigh 



1999, 



2007 



Rozanov et al. 



20061). 



The role of the solar UV radiation for the Earth's climate is boosted 
further by the fact that variations of solar irradiance are also strongly wave- 
length dependent. Whereas the total (integrated over all wavelengths) solar 

ar cycle, th e UV i rradiance 



Floyd et al. 



20031 ). 



irradiance varies by only about 0.1% over the so , 
varies by 1 to 3 orders of magnitude more (e.g.. 

Solar near-IR radiation absorbed by water vapour a nd carbon dio xide is 
an important source of heating in the lower atmosphere (jllaigh . 2007). Solar 
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variability in the IR is comparable to or lower than the TSI variations and 

in the range be tween about 1500 and 2500 nm it is revers ed with respect to 
I II II 



the solar cycle fiHarder et al.l . 12009 : iKrivova et al. 



2009bD. 



Solid assessment of the solar forcing on the Earth's climate is still plagued, 
among other factors, by a shortage of reliable and sufficiently long irradiance 
records. The time series of direct space-borne measurements of solar irra- 
diance is only 3 decades long (in case of spectral irradiance, even shorter). 
Consequently, models are required to extend the irradiance time series to 
earlier times. In the last years, considerable advances have been made in 
modelling the solar total irradiance. Also, spectral models started gathering 
force. Reconstructions of the past variations are quite accurate for the last 
decades, when sufficiently high resolution measurements of the solar pho- 
tospheric magnetic flux (i.e. magnetograms) exist. On longer time scales, 
the secular trend remains rather uncertain and until very recently the recon- 
str uctions were lim i ted to the telescope era, i.e. to the last 4 centuries. 



Domingo et al. 



(120091 ) have recently presented an overview of different 
types of irradiance models on time scales of days up to the solar cycle. Here 
we review the most recent progress in modelling solar total and spectral ir- 
radiance on time scales from days to centuries and even longer using the 
SATIRE set of models. Before discussing the model we briefly describe the 
available observational data (Sect. 12.11) and main sources of irradiance vari- 
ations (Sect. 12. 2p . In Sect. [3] we concentrate on the period when direct 
irradiance measurements are available, in Sects. H] and [5] we discuss the re- 
constructions for the telescope era and the flrst steps towards longer-term 
reconstructions, respectively. Section [6] summarises our current knowledge 
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on irradiance variations on time scales of interest for climate research. 



2. Solar Irradiance Variations 



2.1. Observations 

2.1.1. Total Solar Irradiance 

Measurements of the total solar irradiance (TSI) from space, having suffi- 
cient accura cy to detect i t s vari ations, started with the Hickey-Frieden (HF) 



radiometer (IHickey et al. 



1980) aboard the NOA A/NASA mission Nimbus- 



7. It was followed by t 



le Active Cavity Ra diometer for Irradiance Monitoring 



(ACRIM) I on SMM flWillson et al 



perim ent (ERBE) on ERBS flLee et allllQQSh . ACRIM-II on UARS flWillsonl . 



19811). t he Earth Radiation Bu dget Ex 



994[). Variabih t y of 



rradiance and Gravity Oscillati ons (VIRGO) on SoHO 



dFrohlich et all . Il997h . ACRIM-III on ACRIMS AT dWilsonl. 



Total Solar Irradiance Monitor (TIM) on SORCE (IKopp and Lawrence! . l2005l ) 



2001) and the 



The measurements disclosed that the value previously called the solar 
constant in fact varied at various time scales, from minutes to decades. The 
typical magnitude of variations in the last 3 cycles was about 0.1%, although 
changes up to about 0.3% were recorded on shorter time scales (days), e.g. 
around the Halloween storm at the end of October 2003. 

Since none of the instruments survived over the whole period since 1978 
and each of them suffers from its individual degradation, calibration or other 
problems, a construction of a composite TSI record is quite a challenge. 
Therefore, 3 different composites are now avail able: from the Physikalisch- 
Meteorologisches Observatorium Davos (PMOD; 



Frohlichl . 



20061), the ACRIM 
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team (Willson and Mordvinov. 



Belgique (IRMB 



Dewitte et al 



2003r) and the Institut Royal Meteorologique 



2004h . 



They show certain differences, of which the most important regards the 
presence or absence of an upward secular trend between the minima preceed- 
ing cycles 22 (1986) and 23 (1996). This long-standing argument traces its 
origin back to the unforseen gap between the first two ACRIM experiments 
(July 1989 to October 1991), known as the ACRIM gap. Although 2 other 
experiments, HF and ERBE, continued operati ng during these 2 years, the y 



are claimed to be lower precision experiments by lScafetta and WillsonI (120091 ). 
Moreover, the ERBE made measurements only approximately every second 
week, whereas in September 1989 the output signal of the HF radiometer sat- 
urated and the instrument was switched off for several days. Once back in op- 
eration, it showed an abrupt increase in irradiance by about 0.4 W/ m^. Some- 



what 



19951: 



ater another in c rease of similar magnitude was reported ( iLee et al. 



Chapman et al.l . Il996l ). Depending on whether this change is taken 



into account (as in the PMOD composite) or not (ACRIM), the composite 
time serie s show s either no change in the TSI average level from 1986 to 1996 
( iFrohlichl . l2006l ) or a significant increase (IWillson and Mordvinovl . |2003| ) . 

Another uncertainty of the TSI record regards the absolute level of the 
irradiance. The uncertainty in the absolute values of the TSI is significantly 
higher than that of the relative changes. In particular, there remains an un- 
explained difference of about 4.5 W/m^ betw een SORCE/TIM measurements 



and other currently operating instruments (iKopp et al 



20051 ). A compari- 



son of the TIM measurements to the NIST (National Institute of Standards 
and Technology) reference detectors showed a rather good agreement thus 
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providing support to TIM's lower values. The future Glory mission of NASA 
promises to throw light on this, since the TSI measured with the Glory/TIM 
instrument will, for the first time, be compared directly to the TSI measure- 
ments by the reference cryogenic radiometer under the same vacuum con- 
ditions. For studies of climate change, variations in the TSI are of greater 
importance than the exact absolute value. 

2.1.2. Spectral Irradiance 

The record of spectral irradiance measurements, accurate enough to as- 
sess its variations, goes back to the launch of the Solar Mesospheric Explorer 
(SME) in 1981. Although data between 115 and 302 nm were obtained, their 
accuracy is not sufficiently high due to degradation problems. This led to 
the introduction of the Mg II core-to-wing i ndex, which, being a r a tio of 
the fluxes measur e d at close- by wavelengths (IHeath and Schlesinger 



1986 



Viereck and Pugal . Il999l ) , suffers significantly less from instrument degrada- 



tion. Mg II observations were continued by a number of fo l 



Viereck et al. 



owin g missions 



mm . 



and a composite time series was produced by 

The Upper Atmosphere Research Satellite (UARS) launched in 1991 had 
2 instruments for measurements of solar UV irradiance between approxi- 



mately 1 15 and 420 nm, t 



Brueckner et al. 



(SUSIM; 



ment (SOLSTICE; 



le So lar Ultraviolet Spectral Irradiance Monitor 



19931) and the Solar Stellar Comparison Experi- 



Rottman et al. 



19931). UARS SUSIM and SOLSTICE 



measurements revealed the strong spect ral dependence of the solar irradi- 



Floyd et al. 



20031 ). However, the 



ance variations (see, e.g., the review by 
variability above approximately 300 nm remained rather uncertain due to 
insufficient long-term stability of both instruments. 



Since 1996 measurements have also being taken in three spectral channels 
centere d at 402, 500 and 862 nm with the VIRGO Sun PhotoMeter (SPM) on 



SoHO flFr5hhch et al. 



19971 ). The detectors are, however, not stable enough 



to allow an independent estimate of the comparatively weak longer-term 
trends at these wavelengths. 

Data in a broader spectral range became availab le with th e laun ch of the 
Solar Radiation and Climate Experiment (SORCE; iRottmanl . l2005l ) and the 
Environmental Satellite (ENVISAT) in 2003. SORCE has two instruments 
for measuring solar spectral irradiance: SOLSTICE and the Spectral Irradi- 
ance Monito r (SIM) operating at 115-4 25 nm and 200-2700 nm, respectively 



( ISnow et al. 



2005 



Harder et al.l . 



20051 ). The S Canning Imaging Absorption 



spectroMeter for Atmospheric CHarto graphY (SCIAMAC HY) on ENVISAT 



observes in the range 240-2380 nm (jSkupin et al. 



20051 ) . The data from 



SIM and SCIAMACHY provided, for the first time, information about solar 
variability in the IR part of the sp ectrum and about longer-term (tim e scale s 



of years) variations above 300 nm (IHarder et al. 



2009 



Pagaran et al. 



20091). 



These observations cover only the descending phase of solar cycle 23. The 
data over broad spectral sub-intervals are available, although there are still 
some remaining problems in some others (in particular, in the IR). 

In 2008, the SOLar SPECtrum (SOLSPEC) instr ument on SOLAR/IS S 



20091). 



started observations between 165 nm and 3080 nm (IThuillier et al. 
The data are not yet publicly available. 

Note that here we only discuss the spectral range above 115 nm. Varia- 
tions of the irradiance at yet shorter wavelengths is of more interest for space 
weather studies and their measurements and modelling have recently been 
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reviewed by IWooda (120081 ) . 



2.2. Sources of irradiance variations 

Just a view of the solar surface with a modern telescope (and to some 
extent even without) suggests that brightness of the solar disc is not homoge- 
neous and changes in time. Sunspots are the largest magnetic structures on 
the solar surface and have been seen already in antiquity with the naked eye. 
Being cooler than the surrounding photosphere, they appear darker. Thus a 



passage o: 
(see, e.g.. 



a spot across the visible solar disc caus es a dip in solar brightness 



Hudson et al. 



1982 



Fligge et al. 



2000l ). The spots appear partic- 



ularly dark when they are close to the disc centre when their projected area 
is largest and the dark umbra is best visible. The related irradiance changes 
occur on time scales of hours to about 10 days. 

Since the amount of and area covered by sunspots increases significantly 
from activity minimum to cycle maximum, the frequency and depths of 
brightness dips related to spots also increase considerably. Some of the dips 
reach up to 3-4 W/m^ (or about 0.2-0.3%). On average, however, the Sun 
brightens from cycle minima to maxima. This is because of the contribution 
of faculae and the network , which are weaker magneti c elements appearing 



bright on the solar surface (jZwaan 



1978 



Solanki 



19931 ). Faculae are formed 



in active regions and typically appear in the vicinity of sunspots. The weaker 
network elements are spread over the whole disc, forming a bright network 
at the intersections of supergranules. In the visible continuum, the contrast 
of faculae is significantly higher close to the limb. Thus when a facular re- 
gion crosses the visibl e solar disc i t typic ally causes a weakly double-peaked 



brightening (see, e.g. 



Fligge et al 



2000) 



Wenzlerl ( 120051 ) has considered the net contributions of sunspots and fac- 
ulae to the TSI variation of cycle 23. He found that sunspots caused a 
darkening of about 0.84 W/m^, whereas faculae produced a brightening of 
about 1.73 W/m^. The net increase in the TSI thus was 0.89 W/m^ in 
good agree ment with the value of 0.829±0.017 W/m^ given by the PMOD 
composite (IFrohlichl . |2006| ) . 

Variations in the brightness of the Sun are also produced by the contin- 
uously evolving granulation pattern. Granules are signatures of convective 
cells with sizes ranging from a few hundred to roughly 2000 km. Their cen- 
tral parts, where hot material streams up from the interior to the surface, 
appear brighter, whereas at the outer edges the cooled darker material flows 
down, producing an RMS contrast of appr oximately 14% in the green part of 



the solar spectrum (IDaniloyic et al 



are about 5 to 10 minutes ( 



2008r). The ty pical lifetimes of granules 



Hirzberger et al.l . Il999l ). Granulation dominates 



irradiance variations on time scales shorter t 



lan approxirnately half a day 



and is not important for the Earth's climate (ISolanki et al 



20031) 



2.3. The fundamentals of the SATIRE model 

Here we describe the basic assumptions and ingredients of SATIRE, a 
set of codes to calculate spectral solar irr adiance frorn solar disc images or 



activ i ty proxies such as sunspot n umbers (IFligge et al. 



2003 



1998 



Krivova et al. 



2007 



Wenzler et al. 



20061 ). The assumption underlying SATIRE is 



that all changes in the solar irradiance on time scales longer than hours are 
solely due to changes in the solar surface magnetic flux as traced through 
surface features, such as spots and faculae (see Sec. 12.21 for more detail). To 
obtain the solar irradiance at a given time, t, and wavelength. A, it then 



suffices to derive ffiling factors, /""{t), quantifying how much of the solar 
surface is covered by a given surface feature a (where a denotes umbra, 
penumbra, faculae or network), and to calculate the emergent radiative flux 
of the feature F"(A). The solar irradiance, S{t,X), can then be calculated 
by summing the emergent radiative fluxes of the quiet (inactive) Sun, -^'''(A), 
and those of the surface features according to filling factors 



S{t,X) 



a \ a / 



This formalism can be extended to include positional information, in 
which case the global filling factors f"" are replaced by local filling factors 
a" J, where i,j are indices referring to a particular pixel in a solar image 
or magnetogram. The radiative fluxes are replaced by intensities /"(/x), 
where ^ = cosG and G is the heliocentric angle. Then Eq. ([T]) needs to be 
replaced by 

Sit^ ^) = E E + (1 - ^(..)(^)] • (2) 

Note that the radiative fluxes are independent of time; their calculation is 
discussed in the following paragraphs. The filling factors represent the chang- 
ing aspect of the solar surface, and it is their time dependence that produces 
the variability in the solar irradiance. The filling factors can be obtained 
from a variety of data sets, depending on the time scales considered. The 
type of data available also determines whether Eq. ([2]) is employed (basically 
limited to the satellite era) or Eq. ([T]), for the longer time scales. In order 
to distinguish between models adjusted to different time scales and based 
on different kinds of input data used to derive the filling factors, we call 
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them SATIRE-S (Satellite period, magnetogram-based, Sect. [3]), SATIRE-T 
(Telescope era, from the sunspot number, Sect.H]) and SATIRE-M (Millennia 
time scales, based on cosmogenic isotope data. Sect. [5]). Note, that positional 
information for the photospheric magnetic features is only available in case 
of SATIRE-S, when magnetograms are used to derive the filling factors. 
The emergen t intensi t ies are calcu lated from a set of model atmospheres 



using ATLAS9 (jKuruca . 



1993 



20051 ). This yields relatively low-resolution 



spectra (1 to 2 nm in the NUV and optical), but has the advantage of cov- 
ering the whole wavelength range contributing significantly to the TSI while 
remaining computationally manageable. As ATLAS9 uses an LTE approach 
and combines the opacity contributions due to lines into opacity distribution 
functions, the accuracy of the calculations decreases at shorter wavelength 
and for the strongest lines. For these wavelength regions, more reliable results 
can b e obtained using a correlation analysis as described in iKrivova et al. 

The mo del atmosphe r es ari d resulting emergent intensities are described 



in detail by 



Unruh et al. 



standard solar model of 



1999 ). The m odel for the quiet Sun is based on the 
1 

Kuruca (119921 ): the umbral and penumbral models 
are radiative equilibrium atmosphere models with effective temperatures of 
4500 K and 5400 K, r espectively. Th e facu lar model, finally, is a modified 



version of FAL P from 



Fontenla et al. 



fll993h . 



Solar irradiance reconstructions based on these models have been used 
successfully to calculate the spectral variability over a large range of wave- 
lengths (see Sec. 13. Sp : they also match the measured TSI on solar-cycle 
timescales (Sec. 13. 2p . Despite these successes, a number of issues remain. 
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One of these is that the use of 1-D atmospheres leads to the facular contrast 
being overestimated towards the solar limb. Another is the assumption that 
all small-scale magnetic features are assumed to display the same wavelength 
and limb-dependent contrast behaviour that is scaled according to their mag- 
netic flux. This assumption is unlikely to hold, as suggested by the d i fferent 
thermal structure of network and active region flux tubes ( SolankJ. 1986 ) 



and the magnetic flux d ependence of facular /network contrast (ITopka et al. 



1997; 



Ortiz et al. 



20021 ) . One way to address both of these issues is to use 



realistic simulations of magneto-convection in order to represent the 3-D 
structure of the solar atmosphere and to derive from them the emergent 
intensities as a function of the magnetic flux. 

The flrst steps of such an approach were presented bvlUnruh et aLl (120091) , 



who used MURaM magneto-convection simulations (jVogler et al 



20051) to 



calculate disc-centre intensities. Figure [T] shows flrst results for contrasts at 
/i = 1 (i.e. solar disc centre) derived from magneto-convection simulations 
with mean magnetic fluxes of 50 G, 200 G and 400 G. Also shown is the 
contrast currently used in the SATIRE reconstructions (dotted line); it has 
been reduced by approximately 30% so that it can be compared directly to 
the simulation results. A higher contrast is expected for the SATIRE 1-D 
model as this represents the brightest faculae, while the MURaM contrast 
is averaged over a computational box that contains a mix of magnetic flux 
tubes and fleld-free regions. 

Overall, the spectral behaviour of the 200-G simulations and of the SATIRE 
contrasts is very similar. The 50-G simulations that are more indicative of 
the solar network show a slightly flatter spectral shape with a slower increase 



12 



0.6 h ' ' ' ' ; H 0.03 




P .. I .... I .... I .... I .... ^ -0.04 L ... . . . . . . . j 

250 300 350 400 450 1000 10000 

wavelength [nm] wavelength [nm] 

Figure 1: Disk-centre contrasts of simulated convection with a mean magnetic fiux of 50 G, 
200 G and 400 G compared to a non-magnetic model. The dotted line shows the contrast 
of the 1-D facular atmosphere used in the SATIRE reconstructions. This has been reduced 
by 30% for ease of comparison. 

of the contrast at shorter wavelengths and less darkening in the near-infrared. 



3. Solar Irradiance Over the Satellite Period (SATIRE-S) 



3.1. Temporal Changes from Magnetograms 

The most detailed information on the temporal evolution of the solar 
surface magnetic features, needed in order to calculate irradiance variations, 
is provided by full disc magnetograms. Since 1996 such magnetograms have 
been regularly record ed by the Michelson Doppler Imager (MDI) on SoHO 



( IScherrer et al.l . Il995h . Ground-based magnetograms with a sufficiently large 



spatial resolution, e.g., from the Kitt Peak Nationa . 
KP) are available back to 1974 (ILivingston et al. 



Solar Observatory 



1976 



Jones et al. 



(NSO 



19921), 



although the temporal coverage and the quality become significantly poorer 
for the earlier period. 
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Magnetograms and continuum images recorded as close in time as pos- 
sible are used to identify different magnetic features on the solar surface: 
sunspot umbrae, sunspot penumbrae, faculae and the network, and to deter- 
mine which fraction of the solar surface is occupied by each component at a 
given fi, i.e. the filling factors a (see Sect. 12. 3p . Spot umbrae and penum- 
brae are identified by their continuum contrasts. Pixels belonging to spots 
are then excluded from the (co-aligned) magnetograms and the remaining 
magnetogram signal is considered to be faculae and the network. Pixels with 
the magnetic signal below the noise threshold are counted as quiet Sun. 

The intensity spectra of all components are then weighted by the corre- 
sponding filling factors and summed up to give the irradiance as a function 
of wavelength using Eq. ([2]). Integration over all wavelengths results in total 
solar irradiance at time t. The procedure is repeated for every day when a 
magnetogram and a continuum image are availabl e. The models emplo ying 



MDI and NSO 



Wenzler et al. 



data are descr ibed in detail by 



(12004 



2005 



20061 ). respectively. 



Krivova et a. 



Wenzler et al. 



20m and 



fl2004h also 



compared the reconstructions based on these 2 sets of magnetograms to each 
other. 

The model has a single free parameter, -Bgat? the value of the magne- 
togram signal at which magnetic field in the solar photosphere fills the whole 
magnetogram pixel (i, j). Sgat is found from the best agreement with the TSI 
measurements. Note that one aim of the recently started work on computing 
the spectra from MHD simulations (see Sect. 12. 3p is to do away with this free 
parameter. 
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3.2. Total Solar Irradiance 



Krivova et al.i (i2003ii ) employed MDI magnetograms to reconstruct the 



TSI between 1996 and 2002. Excellent agreement between the PMOD TSI 
composite and t he model was obtained with a cor relation coefficient of Tc = 



0.96 (r2 = 0.92). 



Wenzler et al. 



(12004 



2005 



20061 ) have extended this recon- 



struction back to 1974 using NSO KP magnetograms and continuum images. 
The correlation for the whole period 1974-2003 is somewhat lower than when 
MDI data are used [vc = 0.91, = 0.83), which is not surprising due to the 
significantly lower quality and artefacts seen in ground-based data (especially 
before 1992 when an older instrument was used at the NSO KP) compared to 
those taken from space. Nevertheless, the good agreement between the data 
and the model provides strong support for the model's basic assumption, 
namely that the dominant part of the irradiance variations on time scales 
longer than about a day is caused by the evolution of solar surface magnetic 
features. 



Wenzler et al 



(l2009l ) have also compared the model to all three currently 
available TSI composites: PMOD, ACRIM and IRMB. For this, they have 
varied the value of the single free parameter (-Bsat, see Sect. 13. ip in order 
to find the best fit to each of the composites individually. The overall best 
fit is found with the PMOD composite (r^ = 0.91, slope 0.98), whereas the 
agreement with the other two composites is markedly poorer (rc = 0.79, slope 
0.82 for ACRIM and = 0.82, slope 0.81 for IRMB). This can be judged 
from Fig. [2|, which displays the difference between the TSI reconstructed by 
SA TIREl-S and each of t h e thre e composites. 



Scafetta and WillsonI (120091 ) proposed that these reconstructions can be 
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Figure 2: The top panel shows the reconstructed daily total solar irradiance based on NSO 
KP data for 3528 days between 1978 and 2003, i.e. from the late rising phase of cycle 21 
to the declining phase of cycle 23. The three lower panels show the difference between 
the reconstructed TSI and the PMOD, ACRIM and IRMB composites, respectively. Each 
dot represents a daily value. The horizontal solid lines indicate a difference = 0, the thick 
dashed lines are linear regressions. The dotted vertical lines indicate periods when the 
individual data sets from HF, AI, All & AIII (ACRIM I, II & III), V (VIRGO: Variability 
of solar IRradiance and Gravity Oscillations experiment on SoHO) and D (DIARAD: 
Differential Absolute RADiometer) were used for the composites. From Wenzler et al. 
(2009) 



16 



used to bridge the so-called ACRIM gap (see Sect. I2.1.ip and to create 
a 'mixed' ACRIM-1 - SATIRE - ACRIM-2 composite. They have com- 
pared ACRIM-1 and ACRIM-2 data directly to the model, in order to cross- 
calibrate the data from the two instruments. These authors have, however, 
used the SATIRE-T model described in Sect. IH which is not suited for such 
an analysis. As discussed later, SATIRE-T is based on the historic sunspot 
number record instead of on magnetograms and continuum images, so that 
it is based on Eq. ([T]) rather than the more realistic Eq. ([2]). It is therefore 
significantly less accurate than SATIRE-S on time scales of weeks to sev- 
eral months. These are, however, the most critical time scales for such a 
comparison of the data and the model. This indicates that cautio n need s to 
be exercised when considering the results of IScafetta and WillsonI (120091 ). in 



particular their conclusion about the upward trend between the minima in 
1986 and 1996. 



Krivova et al. 



(j2009af ) have repeated this analysis employing the more ap- 
propriate SATIRE-S model. The 'mixed' ACRIM-1 - SATIRE-S - ACRIM-2 
composite is shown in Fig. [31 It shows no increase in the TSI from 1986 to 
1996, in contrast to the ACRIM composite. Actually, a slight decrease is 
found. The magnitude of this decrease cannot be estimated very accurately 
from such an analysis but it lies between approximately 0.2 and 0.7 W/m^ 
(0.013-0.05%). 



3.3. Spectral Solar Irradiance 

Early comparisons between SATIRE-S spectral irrad iance reconstructions 
and U ARS /SUSIM measureme nts were presented by 



Krivova and Solanki 



([20051) and 



Krivova et al 



(120061 ) for wavelengths below 400 nm. These com- 
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Figure 3: 'Mixed' TSI composite constructed from ARCIM-1 and ACRIM-2 data (black 
dashed line), with the gap bridged using the SATIRE-S (WSK09 ACRIM) model (asterisks 
connected by grey solid line when there are no gaps). The heavy solid line is the 1-year 
smoothed TSI, and the horizontal da shed line shows the level of the minimum preceding 

|2009a|)- 



cycle 22. From 



Krivova et al 
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Figure 4: Normalised variability from May 2004 to July 2004 as measured with 
SORCE/SIM (solid lines) and modelled with SATIRE (dotted hues) for wavelengths span- 
ning 230 nm to 420 nm (left-hand side) and 420 nm to 1600 nm (right-hand side). The 
filled squares indicate the variability measured in the three SoHO /VIRGO SPM bands at 
402 nm, 501 nm and 863 nm. 



parisons showed that there is good agreement between the measured and 
mo delled wavelen g th int egrated flux between 220 and 240 nm. 

(120061 ) then used this wavelength band as a reference to 



Krivova et al. 



derive regression coeflicients measuring the response of those wavelength re- 
gions that are not well represented by the ATLAS9 fluxes. In this way 
SATIRE can be extended down to 115 nm. For wavelengths above 270 nm, 
the sensitivity of UARS/SUSIM is not sufficient to measure solar variability 
reliably, so that modelled fluxes (i.e. those based on model atmospheres) are 
currently more reliable. 

Comparisons in the visible and near-infrared have only become possi- 
ble with the launch of SORCE and EN VIS AT (Sect . [2X^ . First results 



from SORCE/SIM were presented by 



Fontenla et al. 



(120041 ) . Comparisons 



between SIM data and SATIRE-S covering 3 solar rotations in 2004 were 
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Figure 5: Top curves: flux variations for the Mg II band {left-hand side) and for two bands 
centred at 670 and at 1550 nm {right-hand panels). The dotted hues hnking the squares 
and triangles represent SIM data, the plus signs linked by the solid lines show the model 
calculations. Bottom curves: modelled times series for the facular an d spot contributions . 



The upper lines are for the faculae, the lower for the spots. Following lUnruh et al.l (j2008l ). 
carried out by 



Unruh et al 



(120081 ) ■ The measured and modelled variabil- 
ity is shown in Fig. |H There is good agreement between SATIRE-S and 
SORCE/SIM, in particular for wavelengths above 360 nm. Between 310 nm 
and 360 nm detailed comparisons are hampered as the sensitivity of the 
detector is not quite sufficient to measure solar variability adequately and 
the modelled irradiances yield more reliable results. For wavelengths above 
550 nm SATIRE-S shows a slightly flatter decrease in variability than the 
SORCE measurements, though detailed light-curve comparisons show excel- 
lent agreement in the short-term variability. 

Figure [5] shows measured and modelled time series for the Mg II band 
(left-hand side) and for two bands centred at 670 and at 1550 nm, respec- 
tively (right-hand panels). The upper curves in the plots show the measured 
and modelled solar variability while the bottom curves show the modelled 
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contributions of the faculae and spots to the irradiance variations. The fig- 
ures illustrate the change in the lightcurve aspects as one moves from the UV 
to the near-IR. In the UV, as illustrated by the Mg II (at around 280 nm) 
band, the infiuence of the spots is so small that their darkening effect is more 
than compensated for by the faculae, even on solar rotational time scales. 
Furthermore, the facular contrast increase towards the limb is not sufficient 
to counteract the projection effects and the Sun appears brightest when the 
spot groups are at disc centre. At longer wavelengths (e.g., at 670 nm), 
the facular brightenings are much weaker and show a double-peaked struc- 
ture. Thus faculae produce most of the brightening when near the limb and 
very little, if any, when at disc centre. In the near-IR (as illustrated by 
the 1550 nm band), only very little facular brightening is seen close to the 
limb, which is compensated by the earlier onset and longer duration of spot 
da rkening. The overall e ffect of faculae at 1550 nm is one of darkening. 



Krivova et al. 



(j2009bl ) have reconstructed solar UV irradiance in the range 
115-400 nm over the period 1974- 2007 by making use of the empirical ex- 



tension of the SATIRE models by 



Krivova et al 



(120061) employing SUSIM 



data. The evolution of the solar photospheric magnetic flux in this case was 
described by the magnetograms from the KP NSO between 1974 and 2003 
and by the MDI instrument on SoHO since 1996. The gaps in daily data 
were fllled using the Mg II core-to- wing ratio and the solar F10.7 cm radio 
flux. 

Figure [6] shows the reconstructed solar Ly-a irradiance (light grey crosses; 
red in online version), SUSIM measurements (grey dashed line / green) and 



the composite time series (black dotted line / blue) compiled by 



Woods et al 
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Figure 6: Solar Ly-a irradiance since 1974: reconstructed by SATIRE-S (light grey crosses; 
in online version red), measured by the SUSIM instrument (grey dashed line; in online 
version green) and compiled by Woods et al. (2000; black dotted line; in online version 
blue). From Krivova et al. (2009). 
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(120001 ). The latter record includes measurements from the Atmospheric Ex- 
plorer E (AE-E, 1977-1980), the Solar Mesosphere Explorer (SME, 1981- 
1989), UARS SOLSTICE (1991-2001) and the Solar EUV Experiment (SEE) 
on the TIMED mission launched in 2001. The gaps are filled in using proxy 
models based on Mg II and F10.7 indices. The F10.7-based model is also 
used for the period before 1977. All measurements and models are adjusted 
to the SOLSTICE absolute values. The model agrees well with the SUSIM 
data, which confirms that the semi-empirical technique works well. There is 
some difference in the magnitude of the Lj-a solar cycle variation between 
SOLSTICE and SUSIM, which remains also if the model is compared to the 
composite. This is because the model is adjusted to SUSIM by construction. 
Other than that the model agrees with the composite record very well with 
a correlation coefficient of 0.95. 



Employing the solar F10.7 cm radio fiux data (ITanaka et al.l . Il973l ). it 
is possible to extend the reconstruction back to 1947. For this, the radio 
fiux measurements are used to calculate the irradiance at 220-240 nm, in 



the same way as done by 



Krivova et al 



(l2009bl ). who employed these data 



in order to fill in the gaps in their reconstructions for the period 1974-2007. 
Irradiance at other wavelengths is then calculated from the irradiance in 
the re ference range, 220-240 nm, by applying the technique of iKrivova et al. 



(120061 ) described above. Figure [7] shows the reconstructed solar irradiance 
integrated over spectral ranges of particular interest for climate studies as a 
function of time: (a) 130-175 nm, (b) 175-200 nm, (c) 200-242 nm, (d) 200- 
350 nm, (e) 800-1500 nm and (f ) 1500-2500 nm. The model is based on radio 
fiux data before 1974 and on magnetograms after 1974. Solar UV irradiance 
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Figure 7: Reconstructed solar irradiance in the period 1947-2007 integrated over the 
wavelength ranges: (a) 130-175 nm (Schumann-Runge continuum), (b) 175-200 nm 
(Schumann-Runge bands), (c) 200-242 nm (Herzberg continuum), (d) 200-350 nm 
(Hartley-Huggins bands), (e) 800-1500 nm (H2O bands) and (f) 1500-2500 nm (H2O 
bands). 



24 



(panels a-d) varied by a factor 10 to 100 more than the TSI. Interestingly, 
the irradiance in the range between around 1500 nm and 2500 nm (panel f) 
shows an inverse solar cycle variation compared to the TSI and other sp ectral 



Harder et al. 



20091 ) from 



ranges. This inverse variation was also noticed by 
the analysis of SORCE data. Th is is due to the low or even negative contrast 

see also Fig. [5]), so that 



of faculae at these wavelengths (lUnruh et al. 



2008 



their brightening (if any) no longer compensates darkening due to sunspots. 



4. Solar Irradiance Since the Maunder Minimum (SATIRE-T) 

4.I. Temporal Changes from the Sunspot Number 

The advantage offered by high resolution magnetograms of having infor- 
mation on the spatial distribution of the solar surface magnetic features as 
well as the magnetic field strength is not available prior to 1974. Models going 
further back in time have to content themselves with sparser data. One such 
historic record is the sunspot number. It is widely used for reconstructions 
of solar irradiance on time scales longer than a few decades and is, in fact, 
the only direct proxy of solar magnetic activity going back to the Maunder 
mi nimum (a peri o d of unusua lly low solar activity in the 17th century). 



Solanki et al. 



(2000 



2002) have constructed a physical model which al- 



lows a reconstruction of the evolution of the solar photospheric magnetic field 
from the sunspot number. In this way, it is possible to calculate not just the 
total magnetic flux emerging on the solar surface at a given time, but also in- 
dividual contributions of different components: active regions (sunspots and 
faculae), ephemeral regions (smaller short-lived bipolar magnetic regions) 
and the open flux (a small part of the photospheric magnetic flux whose 
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field lines reach into the heliosphere). The calculated open magnetic flux 
is found to be in good agreement with the empirical rec o nstruction of the 



heliospheric magnetic field since 1868 by ILockwood et al. 



(11999 



20091 ) from 



the geomagnetic aa-index. Also, the modelled total photospheric magnetic 



field reproduces the measurements since 19 74 carried out at 



Kitt Peak and Wilcox Solar Observatories (lArge et al. 



: he Mt. Wilson , 



2002 



Wenzler et al. 



2om 



Krivova et al. 



fl2007h and 



Balmaceda et al. 



(120071 ) made use of this model 



in order to also reconstruct solar total irradiance since 1700 and 1610, re- 
spectively. In order to calculate the filling factors (Sect. 12. 3p of sunspots as 
a function of time, t hey employed the sunsp ot areas. These have been mea- 



sured since 1874 (see 



Balmaceda et al 



20091 ) and were derived for the earlier 



period using a linear relationship between sunspot areas and numbers. The 
sunspot contribution is then subtracted from the active region magnetic flux 
to assess the flux present in the form of faculae. The sum of the ephemeral 
region and open magnetic flux describes the evolution of the network. 

The network and facular magnetic fluxes are then converted into corre- 
sponding filling factors (see Sect. 12.31) assuming that the filling factors are 
proportional to the corresponding magnetic field strengths until the satura- 
tion limit {Bsat, Sect. 13.11) is reached. Finally, the irradiance is calculated by 
summing up disc-integrated brightnesses (i.e. radiative fluxes) of different 
components (the quiet Sun, sunspot umbrae, sunspot penumbrae, faculae 
and the network) weighted by their filling factors a s described in Se c t. 12.31 
In addition to the heliospheric flux reconstruction by 



Lockwood et al 



(119991 ) 



and the total magnetic flux data, the model also reproduces well the mea- 
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Figure 8: 11-yr running mean of the reconstructed TSI based on the Zurich (solid line) 
and Group sunspot number (dashed line). 



sured TSI since 1978, although t he accuracy o: 
scales of weeks to about a year (iKrivova et al. 



the model is reduced on time 



2009al . see Sect.ES]). 



4-2. Irradiance 



Krivova et al 



(120071 ) have employed the Zurich sunspot number in order 



Balmaceda et al- 



to rec onstruct solar total irradiance since 1700, whereas 
( 120071 ) used the Group sunspot number in order to go back to 1610. The two 
reconstructions are shown in Fig. [HI They show slightly different irradiance 
levels during the 18th and 19th centuries, but both suggest that, averaged 
over 11 years, the TSI has increased by about 1.3 Wm~^ since the end of 
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the Maunder minimum. By applying different assumptions regardin K the 



Krivova et al 



(120071 ) have 



evolution of the ephemeral region magnetic flux 
also estimated the uncertainty range as 0.9 to 1.5 Wm~^. The lower limit 
i s in g ood agreement with the results of the modelling effort by 



Wang et al 



(120051 ) ■ who neglected the extended length of the ephemeral region cycle, 
which is the main source of the secular trend in the SATIRE-T model. The 
.5 Wm~^ is close to the value of 1.7 Wm~^ from . 



up per limit of 



cf. 



Posted (120041 ) ■ 



LockwoodI (120051 ). who assumed that the network disappeared completely 



during the Maunder minimum. A comparably lo w value for th e magn itude 



Crouch et al. 



(l2008h from 



of the secular increase has also been obtained by 
their model of TSI variations. The model does not, however, go back to the 
Maunder minimum, but stops at 1874. 



All recent estimates of the magnitude of the secular trend are thus signifi 



cant 



1990 



y lower than the va 



White et al. 



1992 



ues deducec 



Lean et al. 



earlier from stellar data (IBaliunas and Jastrow 



19951 ). Nevertheless, even though all re- 
cent models agree that TSI increased by significantly less than 2 Wm~^ since 
the Maunder minimum, considerable uncertainty remains. This is because all 
models have free parameters, or make assumptions that cannot be rigorously 
tested. 

A potential independent source of information on the magnitude of the 
secular trend are the historic archives of full-disc solar images in the Ca II K 
line. Such observations were carried out by a number of observatories around 
the globe since the beginning of the 20th century. Some of the archives have 



recentl y been digitised. Ear l ier analysis o: 



by, e.g., 



Karivappa and Papl (ll996h : lFoukall (119961 ): 



the d i gitised data was pre s ented 



Caccin et al. 



(Il998h : 



Worden et al 
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(Il998l ). Usage of these data is until now limited because of the numerous 
problems and artifacts affecting the images, such as plate defects and aging, 
degradation of instrumentation and obser ving programmes, g e ometrical dis- 



ErmoUi et al 



tortio us, photometrical uncertainties etc. (lErmoUi et all 120071 ) . 
(I2OO9I ) have analysed the quality of the data of 3 digitised historic archives 
from Arcetri, Mt. Wilson and Kodaikanal observatories. They have also been 
able to fix some of the problems, such as, e.g., disc eccentricity. Some others 
are more difficult to correct. These include photometric uncertainties (due 
to, e.g., stray light or missing photographic calibration) and changes in the 
image content with time (due to, e.g., varying wavelength of the observa- 
tions)^ 



ErmoUi et al. 



(120091 ) have developed a method of image calibration in- 
dependent of the calibration exposures and calculated plage coverage as a 
function of time for the 3 analysed archives (Fig. [9]). The Pearson corre- 
lation coefficients for each pair of the data sets lie between 0.85 and 0.93. 
The most significant differences are found for cycles 15, 17 and 19, which 
are partly explained by the stray-light degradations at the beginning of each 
observing programme and by di fferences in the wav elength of observations. 

(I2OO9I ). who also analysed and 



Similar results were obtained by 



Tlatov et al. 



compared Kodaikanal and Mt. Wilson observations (see their Fig. 4). These 
results suggest that caution is needed when using such data without careful 
analysis of their problems and intrinsic instrumental variations. 
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Figure 9: Temporal variation of yearly median values of the plage coverage measured from 
the Arcetri, Kodaikanal, and Mt. Wilson series. The error bars represent the standard 
deviation of measured values over the annual interval; for clarity, they are only shown for 



the M W series. Cycle numbers are given at the top of each cycle. From 



Ermolli et al 



30 



5. Solar Activity Over the Holocene 



Before 1610, no direct proxies of solar magnetic activity are available. The 
indirect proxies at hand are the concentrations of the cosmogenic isotopes 
^°Be and ^^C in natural terrestrial archives. These isotopes are produced 
when high-energy cosmic rays enter the Earth's atmosphere and react with 
nitrogen and other atoms. The flux of the cosmic rays at the Earth is modu- 
lated by the solar open magnetic field, which makes it possible to reconstruct 
the open magnetic flux using the production rates of the cosmogenic radioiso- 
topes. 



Solanki et al. 



(|200J)and 



Usoskin et al. 



( 2004 



20061 ) have combined physics- 



based models for each of the processes connecting ^''C and ^°Be concentration 
with sunspot number and thus were able to reconstructed the sunspot num- 
ber over the whole Holocene. This reconstruction has shown that the Sun is 
currently in a state of unusually high magnetic activity and that prolonged 
periods of very low solar activity similar to the Mau nder minimum (grand 



minima) have regula rly occurred in the past (see also 



Schove 



1955 



1980 



Siscod . 



Krivskyl . Il984l ). Also periods of higher average d solar activity , simi 



Usoskin et al. 



(12007h 



lar to the recent one (grand maxima) have occurred, 
give a list of all grand minima and maxima since 9500 BC detected in the 
reconstructed sunspot number record. 

The reconstruction of solar irradiance for the Holocene is complicated 
by the fact that only a cycle-averaged sunspot number can be derived from 
the ^^C and ^'^Be concentrations. Thus the SATIRE-T model described in 
Sect, m cannot be applied directly. Therefore, first relationships are derived, 
using the SATIRE-T results, between the cycle averaged open magnetic flux 
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(the primary solar quantity obtained from ^^C and ^°Be concentrations) and 
the total magnetic flux, the sunspot number and solar irradiance. These 
relationships are then used to calculate the solar irradiance over the Holocene 
(Vieira et al., in prep.). 



6. Summary 

As weak as they are (about 0.1% over the solar cycle), solar irradiance 
variations affect the Ea rth's climate at a level of 0.2°K over the solar cycle 



(jCamp and Tuna . 120071 ). The role of the total solar irradiance might be sig- 
nificantly reinforced by the wavelength dependence of irradiance variations. 
The amplitude of the variations is 1 to 3 orders of magnitude higher in the 
UV part of the spectrum than in the visible or JR. At the same time, solar 
UV radiation is an important driver of chemical and physical processes in 
the Earth's upper atmosphere. 

Reliable assessment of solar forcing on the Earth's climate is still plagued, 
among other factors, by a shortage of reliable and sufficiently long irradiance 
records. The time series of direct space-borne measurements of solar irra- 
diance is only 3 decades long (in case of spectral irradiance, even shorter). 
Therefore we have recently focussed our main attention on tests of the re- 
constructed irradiance against the newest spectral data, improvement of the 
models in the UV and appraisal of the magnitude of the secular trend in 
irradiance. 

We have reconstructed solar spectral irradiance at 115-160000 nm back 
to 1947. Reconstruction of solar total irradiance goes back to 1610 and 
suggests a value of about 1-1.5 Wm~^ for an increase in the cycle-averaged 
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TSI between the end of the Maunder minimum and the last 50 years, which 
is significantly lower than previously accepted but agrees with other recent 

estimates. First steps have also been made towards reconstructions of solar 
total and spectral irradiance on time scales of millennia. 
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